Competing Phases, Strong Electron-Phonon Interaction and Superconductivity 
in Elemental Calcium under High Pressure 



Z. P. Yin 1 , F. Gygi 2 , and W. E. Pickett 1 
1 Department of Physics, University of California Davis, Davis, CA 95616 
2 Department of Applied Science, University of California Davis, Davis, CA 95616 

(Dated: October 30, 2009) 

The observed "simple cubic" (sc) phase of elemental Ca at room temperature in the 32-109 GPa 
range is, from linear response calculations, dynamically unstable. By comparing first principle 
calculations of the enthalpy for five sc-related (non-close-packed) structures, we find that all five 
structures compete energetically at room temperature in the 40-90 GPa range, and three do so in 
the f00-f30 GPa range. Some competing structures below 90 GPa are dynamically stable, i.e., no 
imaginary frequency, suggesting that these sc-derived short-range-order local structures exist locally 
and can account for the observed (average) "sc" diffraction pattern. In the dynamically stable phases 
below 90 GPa, some low frequency phonon modes are present, contributing to strong electron- phonon 
(EP) coupling as well as arising from the strong coupling. Linear response calculations for two of 
the structures over 120 GPa lead to critical temperatures in the 20-25 K range as is observed, and 
do so without unusually soft modes. 



I. INTRODUCTION 

One of the most unanticipated developments in 
superconducting critical temperatures (T c ) in the 
past few years has been achievement of much higher 
values of T c in elemental superconductors by the 
application of high pressure, and that these im- 
pressive superconducting states evolve from sim- 
ple metals (not transition metals) that are non- 
superconducting at ambient pressure. The first 
breakthrough arose in Li, with T c approaching^ 20 
K, followed by yttrium" a t megabar pressure also 
superconducting up to 20 K and showing no sign 
of leveling off. Both of these metals have electron- 
phonon (EP) coupled pairing, according to several 
linear response calculations^^ of the phonon spec- 
trum, electron-phonon coupling (EPC) strength, 
and application of Eliashberg theory. These impres- 
sive superconductors have been surpassed by Ca, 
with T c as high as 25 K reported^ near 160 GPa. 
Perhaps more unusual is the report, from room tem- 
perature x-ray diffraction (XRD), of a simple cubic 
(hence far from close-packed) structure over a vol- 
ume reduction of 45->30% (32-109 GPa). Whether 
these two unique phenomena are connected, and in 
what way, raises fundamental new issues in an area 
long thought to be well understood. 

Face-centered cubic (fee, Ca-I) at ambient pres- 
sure, calcium transforms at room temperature to 
body-centered cubic (bec, Ca-II) ali^ 20 GPa, is 
identified as simple cubic (sc, Ca-III) in the very 
wide 32-109 GPa range as mentioned above, and 
shows additional phases (Ca-IV, Ca-V) at even 
higher pressures. A sc structure for an element is 



rare, occurring at ambient pressure only in polo- 
nium and under pressure only in a handful of ele- 
mental metalsi 11 ' 12 This identification of a sc struc- 
ture for Ca is particularly problematic, since it has 
been shown by linear response calculations of the 
phonon spectrum by a few group o 13 i 14 ' 15 that (at 
least at zero temperature) sc Ca is highly unstable 
dynamically at all volumes (pressures) in the region 
of interest. Since these calculations are reliable for 
such metals, there are basic questions about the "sc" 
structure itself. 



II. COMPARISON TO RELATED METALS 

Strontium, which is isovalent with Ca, like Ca su- 
perconducts under pressure and undergoes a series 
of structural transitions from close-packed structure 
to non-close-packed structure at high pressure. Sr 
transforms from a fee phase to a bec phase at 3.5 
GPa and then transforms to Sr-III at 24 GPa, to Sr- 
IV at 35 GPa and to Sr-V at 46 GPa^ The Sr-III 
structure was first believed to be a distorted sc and 
later found to be an orthorhombic structure ¥ How- 
ever, later experiments have found that there are 
two phases coexisting in the Sr-III phase, namely, 
a tetragonal phase with a distorted /3-tin structure 
and an unidentified additional phased The Sr-IV 
structure is very complex and was shown recently to 
be a monoclinic structure with the la space group 
and 12 atoms per unit cell^ The structure is more 
complex in Sr-V, and was identified as an incom- 
mensurate structure similar to that of Ba-IV^ Sr 
begins to superconduct at 20 GPa, its T c is 8 K at 58 



GPa, and is believed to be higher beyond 58 GPa. 16 
Scandium, with one more (3d) electron than Ca, 
undergoes phase transitions from hep to Sc-II at 20 
GPa and to a Sc-III phase at 107 GPa^Sl Although 
Sc is conventionally grouped together with Y and 
the lanthanide metals as the rare-earth metals, due 
to their similarities in their outer electron configu- 
rations, its structural transition sequence is rather 
different from the common sequence of lanthanide 
metals and Y, which follow the pattern hep— > Sm- 
type— ► dhep^ fec^ distorted fee. The Sc-II struc- 
ture is complex, and was recently found to be best 
fitted to a pseudo bee structure with 24 atoms in the 
unit celli^ The structure of Sc-III is not identified 
to date. Sc begins to superconduct at 20 GPa. Its 
T c increases monotonically to 19.6 K with pressure 
to 107 GPa. Its T c drops dramatically to 8 K at 
the phase transition from Sc-II to Sc-III around 107 
GPa^i 

Considering the close relation of Sc and Sr to Ca 
in the periodic table and the similar superconduct- 
ing properties under pressure, it could be expected 
that Ca under pressure should have more complex 
structures, rather than the observed sc structure. In 
fact, Olijnyk and HolzapfeU^ observed that their Ca 
sample transformed from sc to an unidentified com- 
plex structure at 42 GPa. 

So far the higher pressure phases Ca-IV and Ca- 

V have attracted the most attention, and consid- 
erable progress has been made in identifying these 
phases through a combination of experimental ^ 22 ' 23 
and theoretica l 24 ' 25 ' 26 work. However, satisfac- 
tory agreement between experimental and theoret- 
ical work is still lacking. Ca-IV is identified as a 
Pnma space group by Yao et al^ but P4a2i2 sym- 
metry by Ishikawa et al£2. and Fujihisa et alM- Ca- 

V seems clearly to have a Cmca space group 23 ' 24 ' 25 , 
however, the calculated enthalpy in the Pnma struc- 
ture is much lower than in other structures (includ- 
ing Cmca structure) at pressures over 140 GPa. Also 
in the experimental work of Fujihisa et o/i^ 3 -, the fit- 
ting of their XRD patterns to the anticipated P4a2i2 
and Cmca space groups were not satisfactory and 
other possibilities still exist. In the recent work 
of Arapan, Mao, and Ahuja— , an incommensurate 
structure similar to Sr-V and Ba-IV structures was 
proposed for Ca-V phase. Therefore the nature of 
the Ca-IV and Ca-V phases is still not fully settled. 

While helping to forge an understanding the struc- 
ture of Ca-IV and Ca-V and its impressive supercon- 
ducting T c is one goal of the present work, our focus 
has been to understand the enigmatic "sc" Ca-III 
phase where relatively high T c emerges and increases 



TABLE I: Detailed structural data of the I43m, Pnma, 
Cmca and P4 3 2i2 Ca. (SG: space group; WP: Wyckoff 
position; AC: atomic coordinates.) 



SG 


No. 


WP 


AC 


X 


y 


z 


I43m 


217 


8c 


i^x 2 X j X^j 


~ 0.2 






Pnma 


62 


4c 


(x, 1/4, z) 


~ 0.3 




~ 0.6 


Cmca 


64 


8f 


(0, y, z) 




~ 0.3 


~ 0.2 


P4.-A2 


96 


8b 


O, y, z) 


~ 


~ 0.3 


~ 0.3 



with pressure, a phase that XRD at room temper- 
ature (Tjj) identifies as primitive simple cubic. 22 In 
this pressure range sc Ca becomes favored over the 
more closely packed fee and bec structures, but the 
dynamical (in)stability was not calculated by Ahuja 
et aliU We report here first principles calculations 
of the enthalpy of five crystal structures (with space 
groups sc, I43m, P432i2, Cmca, and Pnma), and lin- 
ear response calculations of EPC, that helps to clar- 
ify both the structural and superconducting ques- 
tions. 



III. THEORETICAL APPROACH 

A. Competing Structures 

The most unstable modes of sc Ca are transverse 
[001]-polarized zone boundary modes along the (110) 
directions. A linear combination of the eigenvectors 
of this mode at different zone boundary points leads 
to a body-centered four-atom cell in the space group 
I43m, whose local coordination is shown in the cubic 
cell in the inset of Fig [TJ and has a clear interpre- 
tation as a buckled sc lattice. This structure, when 
relaxed, has no dynamical instabilities. 

The I43m structure is just one kind of distor- 
tion from the sc structure. There are many kinds 
of other possible distortions. Actually several other 
structures including Pnma, Cmca and P4a2i2 were 
proposed for the high pressure Ca-IV and Ca-V 
phases . 23 ' 24 ' 25 ' 26 Their structural details are listed in 
Table IT] and their structures are pictured in reference 
[23ll24ll25l | . I43m is a body-centered cubic structure, 
Pnma and Cmca Ca are orthorhombic, and P4a2i2 
has a tetragonal symmetry. All are closely related 
to sc structure. For example, I43m turns to simple 
cubic if a;=0.25, and the Cmca structure becomes a 
sc structure if a = b = c and y = z = 0.25. 



B. Calculational Methods 

We have used the full potential local-orbital 
(FPLO) code^ the full-potential linearized aug- 
mented plane-wave (FPLAPW) + local orbitals 
(lo) method as implemented in WIEN2K, 29 , the 
Qbox code^S and the PWscf code 3 -^ to do vari- 
ous structural optimizations and electronic struc- 
ture calculations, and check for consistency among 
the results. For the enthalpy calculations we used 
the PWscf code.— Both Qbox and PWscf use 
norm-conserving pscudopotcntials, while the FPLO 
and WIEN2K codes are all-electron and full po- 
tential codes. The linear-response calculations of 
phonon spectra and electron-phonon spectral func- 
tion a 2 F{oj) were done using the all-electron, full 
potential LMTART code^^ 

The parameters used in PWscf for the structural 
optimizations and enthalpy calculations were: wave- 
function planewave cutoff energy of 60 Ry, density 
planewave cutoff energy of 360 Ry, k mesh sam- 
plings (respectively, number of irreducible k points) 
24*24*24 (455), 32*32*32 (897), 24*24*8 (455), 
24*24*24 (3614), 24*32*32 (6562) for sc, I43m, 
P432i2, Cmca, and Pnma structure, respectively. 
Increasing the number of k points lowers the en- 
thalpy by only 1-2 meV/ Ca almost uniformly for 
all structures, resulting in negligible change in vol- 
ume, lattice constants, and internal coordinates. In 
these calculations, we used a Vanderbilt ultrasoft 
pseudopotential 34 with Perdew-Burke-Ernzerhof 35 
(PBE) exchange correlation functional and nonlin- 
ear core-correction, which included semicore 3s3p 
states as well as 4s3g? states in valence states. 



IV. THE ENTHALPIES 

We have calculated enthalpy H(P) curves for 
each structure in the pressure range 40-220 GPa 
based on density functional method a 36 ' 37 using the 
PWscf code^i Several energy differences and re- 
laxations were checked with the Qbox^S FPLOj2& 
and WIEN2K22 codes. In the 40-70 GPa range, 
all five of the structures we have studied have en- 
thalpies that differ by less than 20 meV/Ca (230 
K/Ca), as shown in Fig. H In the 80-100 GPa 
range, the P4s2i2 phase is marginally the more sta- 
ble phase. Three phases are degenerate, again within 
20 meV/Ca, in the 100-130 GPa region and are al- 
most exactly degenerate around 110-115 GPa. Thus 
at room temperature all five phases, including the 
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FIG. 1: (color online) Local coordination of the five 
structures of Ca, plotted as number of neighbors ver- 
sus the distance d relative to the cubic lattice constant 
a sc with the same density. The inset shows the unit 
cube of the I43m structure (which contains two primitive 
cells); this structure retains six near neighbors at equal 
distances but three different second neighbor distances. 
The P4s2i2 and Pnma structures can be regarded to be 
seven-coordinated, albeit with one distance that is sub- 
stantially larger than the other six. 



sc one, are thermodynamically accessible up to 80- 
90 GPa, above which the sc and I43m structures 
become inaccessible. The other three phases remain 
thermally accessible to 130 GPa. Above 140 GPa, 
the Pnma phase becomes increasingly more stable 
than the others. 

Our results agrees well with the results reported 
recently by Yao et alw^ and Ishikawa et al^- in 
their corresponding pressure range. At low pres- 
sure, our result is apparently different from the re- 
sult by Arapan, Mao, and Ahuja^ In their results, 
sc Ca has the lowest enthalpy from 40 GPa to 77 
GPa, lower than the P4a2i2 and Cmca structures. 
A possible reason is that the authors might not have 
taken into account the change in shape and inter- 
nal coordinates of the Cmca structure in the 70-80 
GPa pressure range. In our calculation, b/a=1.0003 
and internal coordinates y=0.254, z=0.225 at 70 
GPa (and similarly below) change dramatically to 
b/a=1.0594, y=0.349 and z=0.199 at 80 GPa (and 
similarly above). 

Although equally dense, quasi-degenerate, and re- 
lated to the sc structure, these structures differ 
in important ways from the sc structure and each 
other. In Fig. Q] the distribution of (first and sec- 
ond) neighbor distances d, relative to the sc lattice 
constant a sc , are pictured. The collection of dis- 
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FIG. 2: (color online) Plot of the enthalpy H(P) of the 
four distorted Ca structures relative to that for Ca in 
the simple cubic structure. The inset gives an expanded 
picture of the 40-100 GPa regime. 



tances cluster around d/a sc ~ 0.97— 1.05 and, more 
broadly, around y/2. In an ensemble of nanocrys- 
tallites of these phases, the radial distribution func- 
tion in the simplest picture should look like a broad- 
ened version of the sc one. For Ca the actual mi- 
croscopic configuration at room temperature, where 
fluctuations (spacial and temporal) can occur among 
these phases (whose enthalpies differ by less than 
&bTr per atom), will no doubt be much more com- 
plex. However, this simplistic radial distribution 
plot makes it plausible that the resulting thermal 
and spatial distribution of Ca atoms will produce an 
XRD pattern more like simple cubic than any other 
simple possibility. Teweldeberhan and Bonev have 
noted the near degeneracy of some of these phases 
in the 40-80 GPa region, and suggest that the T=0 
structure is Pnma in the 45-90 GPa ranged, which 
is consistent with our results if the P4 3 2i2 structure 
is not included. 



V. STABILITY AND LATTICE DYNAMICS 

The structural stability of the (quasi-degenerate) 
structures we have studied provide insight into be- 
havior of Ca under pressure. Linear response calcu- 
lations were performed using the LMTART cod o 32 ' 33 
to evaluate EPC. 

60-100 GPa. The I43m and Pnma structures are 
mostly dynamically stable from 60-100 GPa accord- 
ing to our linear response calculations, but there are 
very soft zone boundary modes that verge on insta- 
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bility (small imaginary frequencies) at some pres- 
sures. The Cmca and P4a2i2 structures are unstable 
over this entire pressure range; note that their struc- 
tures are close to the sc structure. However, they are 
close to stable with very soft phonons at 100 GPa, 
where they were distorted far enough from the sc 
structure. 

A rather common feature among these structures 
in this pressure range is softening of modes at the 
zone boundary, with associated low frequency weight 
in the spectral function a 2 (uj) that can be seen in 
Fig. [3] and Fig. |4j Such low frequency weight con- 
tributes strongly to A, though the contribution to 
T c is better judged^ by < lu > A or even < to 2 > A. 
With increase of pressure, the peaks move towards 
lower frequency, A increases, and the structures ap- 
proach instability. These results are consistent with 
the changes of structure parameters we obtain in the 
process of calculating the enthalpies, where all four 
structures evolve further from the sc structure with 
increase of pressure. 

Above 100 GPa. At the highest pressures studied 
(by us, and experimentally), the crystal structures 
deviate more strongly from the sc structure. Of the 
structures we have considered, the P4a2i2 one be- 
comes favored and also is structurally stable around 
110 GPa. This stability is consistent with the ob- 
served transition from the sc structure to the Ca-IV 
structure at room temperature. The dramatic drop 
in the electrical resistance at around 109 GPa is also 
consistent with a transition from a locally disordered 
phase to a crystalline material^ 

In the pressure range of 110-140 GPa, the 
P432i2, Cmca, and Pnma structures become quasi- 
degenerate again. Linear response calculations of 
the Pnma structure at 120 GPa and above, and of 
the Cmca structure at around 130 GPa indeed show 
strong coupling with A > 1.0 in all the cases. Unlike 
what was found below 100 GPa, there are no longer 
very low frequency phonons (see Fig. [5] and Fig. 
[5]). The coupling strength is spread over frequency, 
peaking for mid-range frequency phonons. 

Another interesting feature arises in the a 2 (us) 
curves, which reveal that the coupling matrix ele- 
ments become relatively uniform across most of the 
frequency range (except the uninteresting acoustic 
modes below 2 THz) at pressures over 120 GPa in 
Pnma structure and at 130 GPa in Cmca structure; 
this behavior is evident in Fig. 2] and especially in 
Fig. O where the results for the Cmca structure at 
130 GPa are pictured. This characteristic is fun- 
damentally different from that below 100 GPa, dis- 
cussed above. 
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At pressures over 140 GPa, the Prima structure 
is clearly favored in our calculation, and linear re- 
sponse calculations indicate the structure is dynam- 
ically stable. The overall results are evident in 
Fig. [4l which shows that the structures remain sta- 
ble (no imaginary frequencies) and the lattice stiff- 
ens smoothly with increasing pressure, and in Fig. 
[5] that shows that strong electron-phonon coupling 
persists and T c remains high. In this high pressure 
range, the incommensurate structure proposed by 
Arapan, Mao, and Ahuja 2 ^ at pressure over 130 GPa 
is also a possibility. 
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FIG. 3: (color online) Plot of a 2 F(u>) (lower panel), 
a 2 {uj) (middle panel), and phonon DOS (upper panel) 
of I43m structure at about 61, 71, 83 and 97 GPa. This 
regime is characterized by strong coupling a 2 (io) at very 
low frequency. 



VI. COUPLING STRENGTH AND T c 



Figure [5] shows the calculated A, r] = Mc a < u 2 > 
A, and rms frequency < uj 2 > x / 2 versus pressure for a 
few structures and pressures. The calculated values 
of T c are shown in the lower panel, using two values 
of Coulomb pseudopotential /i*=0.10 and 0.15 that 
bracket the commonly used values and therefore give 
an indication of the uncertainty due to the lack of 
knowledge of the value of /i* and its pressure depen- 
dence. Results are provided for Ca in I43m, Pnma 
and Cmca structures at a few pressures up to 220 
GPa. In elemental metals and in compounds where 
coupling is dominated by one atom type, r\ has often 

V 



39 



been useful in characterizing contributions to T, 
increases with pressure monotonically by a factor of 
more than 5 from 60 GPa to 220 GPa. The coupling 
constant A increases modestly up to 120 GPa then 
remains nearly constant at A = 1.2-1.4. As pointed 
out elsewhere^ a dense zone sampling is needed 
to calculate A accurately, so any small variation is 
probably not significant. The increase in 77 beyond 
120 GPa correlates well with the lattice stiffening 
(increase in < uj 2 >) in this pressure range. 

The trend of the resulting T c generally follows, 
but seems to overestimate somewhat, the experi- 
mental values 9 . For Cmca structure at about 130 
GPa, the calculated EPC strength is A = 1.2 and 
T c = 20-25K (for the two values of n*) in very sat- 
isfactory agreement with the observed values of T c 
in this pressure range. For Pnma structure, T c in- 
creases rapidly in the 80-120 GPa region. At pres- 
sures above 120 GPa up to the maximum 220 GPa 
that we considered, the EPC constant A is ~ 1.2-1.4 
and the calculated T c increases modestly from 25- 
30K at 120 GPa to 30-35K at 220 GPa. Neither the 
structure dependence nor the pressure dependence 
seems very important: the strong coupling and high 
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FIG. 4: (color online) Plot of a 2 F(Lu) (bottom panel), 
a 2 (uj) (middle panel), and phonon DOS (upper panel) 
of Pnma structure at about 60, 85, 120, 160, and 200 
GPa. The main trends are the stiffening of the modes 
with increasing pressure, and the retention of coupling 
strength a 2 (u>) over a wide frequency range. 



FIG. 5: Phonon spectrum, phonon DOS, a 2 and a 2 F of 
Cmca Ca at 0.251 V (~ 130 GPa from PwSCF). The 
high symmetry points are T(0, 0, 0), Yi(l, 0, 0), Y 2 (0, 
1, 0), F'(l, 1, 0), S(0.5, 0.5, 0), Z(0, 0, 0.5), Ti(l, 0, 
0.5), T 2 (0, 1, 0.5), Z'(l, 1, 0.5) and R(0.5, 0.5, 0.5) in 
the units of (27r/a, 27r/b, 2n/c). 



T c is more the rule than the exception. Ca at high 
pressure may be an excellent superconductor regard- 
less of its structure. 



VII. SUMMARY 

Calculations of enthalpy versus pressure for five 
crystalline phases of Ca (simple cubic and four dis- 
tortions from it) indicate quasi-degeneracy, with en- 
thalpy differences small enough that one might ex- 
pect a locally disordered, highly anharmonic, fluctu- 
ating structure at room temperature. Over most of 
the 30-150 GPa range, we find at least three crys- 
tal phases whose enthalpies indicate they will com- 
pete strongly at room temperature. The sc phase 
itself is badly unstable dynamically (at T=0), but 
the observed "sc" diffraction pattern can be under- 
stood as a locally noncrystalline, highly anharmonic 
phase derived from a spatially inhomogeneous and 
dynamically fluctuating combination of these struc- 
tures, with most of them being straightforward dis- 
tortions from the sc structure. Such a scenario seems 
to account qualitatively for the XRD observations of 
a "sc" structure. 

At pressures below 100 GPa, the quasi-degenerate 
structures tend to have soft branches or occasionally 
lattice instabilities, which are associated with strong 
electron-phonon coupling. In the pressure range of 
110 to 130 GPa three phases (P4 3 2i2, Cmca and 
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FIG. 6: Upper panel: Calculated electron-pfionon cou- 
pling constant A, r\ and To of Ca in I43m (empty sym- 
bols), Pnma (filled symbols) and Cmca (crossing- line 
filled symbols) structures at a few pressures. Lower 
panel: T c calculated from the Allen-Dynes equation, 
showing the dependence on the Coulomb pseudopoten- 
tial for which two values, /i*=0.10 and 0.15 have been 
taken. 
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As our other main result, we find that linear re- 
sponse calculation of the EPC strength and super- 
conducting T c accounts for its impressive supercon- 
ductivity in the high pressure regime and accounts 
in a broad sense for the strong increase of T c in the 
"sc" phase. At higher pressure beyond the current 
experimental limit (i.e., 161 GPa), T c still lies in the 
20-30 K range for some phases that we have studied. 
In fact strong electron-phonon coupling seems to be 
present in several phases across a substantial high 
pressure range, although we have no simple picture 
why such strong coupling should arise. (The strong 
coupling in Li and Y likewise has no simply physical 
explanation.—^) These results may resolve some of 
the perplexing questions on the structure and record 
high T c for an element, and should help in obtaining 
a more complete understanding of the rich phenom- 
ena that arise in simple metals at high pressure. 

After submission of our manuscript, we became 
aware of a study by Yao et al^. They performed 
structural studies of calcium in the range 34-78 GPa 
using mctadynamics and genetic algorithm meth- 
ods. Their methods and results are complementary 
to ours, with each approach providing its own in- 
sights. Connections of their work to ours is evi- 
dent; for example, the I4i/amd structure they fo- 
cused on is slightly distorted from simple cubic, as 
are the structures that we study. Since its enthalpy 
is within 20 meV/Ca of the Pnma structure across 
this pressure range, their result is consistent with 
our explanation of the observation of the simple cu- 
bic diffraction pattern at room temperature. Their 
linear response calculations of electron-phonon cou- 
pling and the resulting T c are also consistent with 
the more extensive results that we present. 
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